The 2-D electron gas at arbitrary spin polarizations and arbitrary coupling strengths: 
exchange-correlation energies, distribution functions and spin-polarized phases. 
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We use a recent approach [Phys. Rev. Letters, 84, 959 (2000)] for including Coulomb interactions 
in quantum systems via a classical mapping of the pair-distribution functions (PDFs) for a study of 
the 2-D electron gas. As in the 3-D case, the "quantum temperature" Tq of a classical 2-D Coulomb 
fluid which has the same correlation energy as the quantum fluid is determined as a function of 
the density parameter rs- Spin-dependent exchange-correlation energies are reported. Comparisons 
of the spin-dependent pair-distributions and other calculated properties with any available 2-D 
quantum Monte Carlo (QMC) results show excellent agreement, strongly favouring more recent 
QMC data. The interesting novel physics brought to light by this study are: (a) the independently 
determined quantum-temperatures for 3-D and 2-D are found to be approximately the same, (i.e, 
universal) function of the classical coupling constant F; (b) the coupling F increases rapidly with 
in 2-D, making it comparatively more coupled than in 3-D; the stronger coupling in 2-D requires 
bridge corrections to the hyper- netted-chain method which is adequate in 3-D; (c) the Helmholtz 
free energy of spin-polarized and unpolarized phases have been calculated. The existence of a spin- 
polarized 2-D liquid near ~ 30, is found to be a marginal possibility. These results pertain to 
clean uniform 2-D electron systems. 

PACS Numbers: 05.30.Fk, 71.10.-hx, 71.45.Gm 



The past three decades have established the two- 
dimensional electron gas (2DEGi as a goldmine of new 
physics and novel technologies.B The 2DEG, ranging 
from a classical to a quantum system, occurs on liquid- 
He surfaces, semiconductor hetero-interfaces, and in the 
CuO layers of high-Tc superconductors. The interest in 
the quantum Hall effect, the metal-ipsulator transition, 
Kosterlitz-Thouless type transitionsjj and applications 
to nanotechnology and spintronics widen its importance. 
The exchange-correlation energy {Exc) of the 2DEG is 
important in its own right as an input to the density- 
functional theory (DFT) of inhomogeneous 2-D systems. 

The physics of the 2DEG depends crucially on the 
"couphng parameter" F = (potential energy) /(kinetic en- 
ergy) arising from the Coulomb interactions. The F for 
the 2DEG at T = and mean density n is equal to the 
mean-disk radius Vg = (tto)"^/^ per electron. In the ab- 
sence of a magnetic field or disorder effects, the density 
parameter r^, the spin polarization C, and the tempera- 
ture T are the only relevant physical variables. A major 
problem in the 2DEG is the need for a theoretical method 
applicable at arbitrary values of F, Q and T . As in the 
3-D electron gas (3DEG), various perturbation methods, 
equatioiaa-of- motion methods, e.g., those of Singwi et al 
(STLS),eI work for relatively weak coupling, but lead to 
negative pair-distribution functions (PDFs) even at mod- 
erate Tg. In spite of a large effort, the calculation of 
the 2DEG gij{r), where i,j are spin species, at arbi- 
trary spin-polarization rg, and T is unsolved, except 
by direct quantum simulation techiuques. Such "quan- 
tum Monte Carlo" (QMC) methods,Q as well as Feenberg 



techniques]^ assume a trial wavefunction ip = FD where 
_F is a correlation factor and _D is a Slater determinant. 
The variational Monte Carlo (VMC) methodp and the 
fixed-node Green-function Monte-Carlo method (GFMC) 
have been applied to the 2DEG, providing the exchange- 
correlation energy, Exc{ts) at T = in a parametrized 
form.Q Finite-T systems involve excited states as well as 
the ground state. While it is easy to get good E^cifs) at 
T — Q — Q, the same is not true for most properties like 
the PDFs, local-field corrections (LFCs), and xc-energies 
at finite T, spin polarization Q and at strong coupling. 

Recently we presented a computationally simple, con- 
ceptually novel classical method for calculating the PDFs 
and other properties (e.g, static response) of the 3DEG al 
arbitrary coupling, spin polarization and temperature.u 
The method is based on identifying a "quantum tempera- 
ture" Tq such that the correlation energy of a correspond- 
ing classical Coulomb fluid at Tq is equal to that of the 
quantum fluid at T = 0. Our only "many-body" input is 
a set of values of Ec{rs) for determining Tq. Once Tq{rs) 
is known, many properties inaccessible to standard meth- 
ods become readily calculable. A discussipja-of the T = 
and flnite-T 3DEG was presented earlier 

A method which works in 3-D does not necessarily 
work in reduced dimensions. A first objective of this let- 
ter is to show that the same ideas apply to the challeng- 
ing problems of the 2DEG. In addition, this study reveals 
previously unexpected aspects of the quantum tempera- 
ture Tq, exposes the long-range nature of the Pauli ex- 
clusion effect in the 2-D system, and provides new results 
for the partially polarized 2-D electron liquid. It is found 



1 



to be numerically precise enough, even at-|rs 

sE 



30, to 

selectively favour newer QMC calculationslH We also ex- 
amine the stability of spin-polarized phases of the 2DEG 
at finite-T and sufficiently large r^. 

The PDFs of the 3DEG could be accurately calculated, 
using the hypcr-netted-chain (HNC) approximation.!!!! 
The 2-D system requires that bridge-type cluster correc- 
tions to the HNC approximation be included. A compar- 
ison of the 2DEG and 3DEG forms of Tq{rs) reveals that 
the quantum temperature is approximately a universal 
function of the classical- fluid coupling constant. 

The essence of our method is to start from the quan- 
tum mechanical g°(r) of the non- interacting problem and 
build up the interacting g[r) by classical methods. 

Consider a fluid of mean density n containing two spin 
species with concentrations Xi = rii/n. We deal with the 
physical temperature T of the 2DEG, while the temper- 
ature Tcf of the classical fluid is 1//3. Since the leading 
dependence of the energy on temperature is quadratic, 
we assume that Tcf= ^/i^"^ + Tg). This is clearly valid 
for T = and for high T, and was justified in more de- 
tail in ref 10, In this letter the main effort is to study 
the 2DEG near T = and at zero magnetic field by de- 
termining Tq. The equations for the PDFs of a classical 
fluid, and the Ornstein-Zernike(OZ) relations are: 

9ij{r) = exp[~l3(j)ij{r) + hij{r) - Cy(r) + By (r)] (1) 

hj{r) ^ Cij{r) + 'Esns J dr'hi^s{\r - r'\)cs^j{r') (2) 

Here <j)ij{r) is the pair potential between the species 
?, J. For two electrons this is just the Coulomb potential 
ycou{f)- If the spins are parallel, the Pauli principle pr&=| 
vents occupation of the same spatial orbital. As before,L! 
we introduce a "Pauli potential", 7^(r). Thus 4>ij{r) be- 
comes P{r)Sij + Vcouir)- The function h{r) — g{r) — 1 is 
related to the structure factor S{k) by a Fourier trans- 
form. The c(r) is the "direct correlation function (DCF)" 
of the OZ equations. The Bij{r) in Eq. 1 is the "bridge" 
term due to certain cluster interactions. If this is ne- 
glected, Eqs. 1-2 form a closed set defining the HNC 
approximation. The HNC is sufficicmt for the 3DEG for 
the range of Vg studied previously.l3 The classical cou- 
pling constant of the 2DEG is found to increase with Ts 
more rapidly than in the 3DEG. Thus, although V{r) 
restricts clustering effects in the case of gn, bridge con- 
tributions (which contain irreducible 3-body and higher 
terms) are found to be important for 1712. The inclusion 
of a bridge term into_Lhe HNC equations, usually via a 
"hard-sphere" model,!!3 is well known in 3-D problems, 
but is less unexplored for 2-D systems. In this context, it 
is noteworthy that Kwon et al. had tg-includc three-body 
correlations in QMC trial functions.EHa 

In the non-interacting system at temperature T, n = 
ni + 712, Xi = rii/n, the anti-|| gi2{r, T) is unity while 
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Here k, r are 2-D vectors and n{k) is the Fermi oc- 
cupation number at the temperature T. At T — Q 
f{r) — 2Ji{kir)/kr where Ji{x) is a Bessel function. The 
Pauli exclusion potential is defined by: 



PV{r)^hUr)~cUr)-ln[gUr)] 



(3) 



where, e.g., c° (r) is the spin-|| DCF of the 0-Z equation. 
At T = it may be shown that (3V{r) ^ —2ln{r) for 

r ^ 0, PVir) - TT/rkp for r -> cx), /?'P(fc) 47r/fc2 

for k — > 00. Only the product pV{r) is determined. The 
classical fluid "temperature" Tcf =1/(3 is still undefined. 
The Pauli potential is a universal function of rkp, and 
limited to about a thermal wavelength Xth at finite T. 

The next step in the method is to use the full (f)ij (r) , 
and solve the coupled HNC and OZ equations for the 
binary (up and down spins) interacting fluid. The 
Coulomb potential Vcoui^) for two point-charge electrons 
is 1/r. However, an electron at the temperature T is lo- 
calized to within a thermal wavelength. ThuSj for the 
3DEG wc used a "diffraction corrected" formrJ, 



Vcouir) = (l/r)[l 



-rkth] 



(4) 



where fct/|_was taken to be fc°^ = (27rm*rc/)^^^ as in Mi- 
noo et al.llj Here m*=l/2 is the reduced mass of the elec- 
tron pair. In the case of the 2DEG we have determined 
kth by numerically solving the Schrodinger equation for 
a pair of 2-D electrons in the potential 1/r and calcH-, 
lating the electron density in each normalized state.o 
Assuming that the diffraction correction has the same 
form as before, i.e., Vd{r)=exp{—kthr) /r, it was found 
that kth/k^^ = 1.158T0-1°3^ where T is in au. 

As in the 3DEG, we need the "quantum temperature" 
Tq for each of the 2DEG, such that the classical 2D- 
fluid has the same Ec{rs) as the quantum fluid. Unlike 
in the 3DEG, an additional complication in the 2DEG 
is the need for a bridge function. In the ||-spin case, 
if the hard-disk radius of the bridge function is smaller 
than the range of the Pauli potential, we may neglect the 
bridge term. In the anti-|j case, we need a bridge term 
B12 explicitly included, say, via a hard-disk model.E£! 
The same bridge term enters the paramagnetic PDF, 
g{r) = 0.5((7ii -|- (712) via the gi2. We can decouple the 
determination of Tq and B12 by obtaining Tq from the 
fully spin polarized system as the bridge function is neg- 
ligible in this case. Using the analytic fit of Ec{rs, C = 1) 
given by Tanatar and Ceperleyjj t = Tq/Tp is found to 
be fitted by: 



t = 2/(1 + G.86413(ry'^ - 1)^) 



(5) 



h^iir) = 2^ki,k2'^(fci)"-(fc2)e 



;(ki 



= -Ifir)? 



This looks very different from the Tq(rs) mapping for the 
3DEG. However, when Tq is plotted as a function of the 
classical coupling constant Tc = l/(Tqrs), (cf. Table I), 
the 3DEG and 2DEG maps are seen to be almost a single 
universal function (Fig. l.(a)). This result is conceptaa/Zj/ 
very interesting since Tq is really a single-parameter rep- 
resentation of the density-functional correlation energy 
of the electron gas, irrespective of dimensionality. 
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Unlike in the 3DEG, we need the bridge function for 
gi2 in order to implement the method. We include only 
the anti-jl Bi2{r), and hence a one-component hard-disk 
packing-fraction 77 is needed. The effective coupling con- 
stant, Teff — (j){rs)/Tcf, is a function of VsTq and ktufs- 
That is, rjt and t^/r^ where t = Tg/Tp. The Gibbs- 
Bogoliubov expression for the free energy of the one- 
component hard-disk system has the same dependence. 
Also, at large the packing fraction should tend to a 
finite limit. Since t decreases as Vs we consider the 

1 /3 

simple functional form rj — Otrg ■ Here 6* is a "global" 
parameter valid for the whole range of and hence the 
procedure does not involve fitting at each Ts ■ A constant 
value, 6 — 0.1175, proved to be applicable to the range 
Ts =1 — 30 studied here. Values of the packing fraction 
77, the QMC and classical-map HNC (CHNC) values of 
Ec(,i"sX = 0) for the unpolarizcd 2DEG are given in Ta- 
ble I. Our Ec{rsX — 1) agrees with Tanatar-Ceperley 
spin-polarized QMC by construction. It is common for 
ry to approach ~ 0.49 at the liquid-solid transition. The 
trend in the value of 77 seen in Table I is consistent with 
the finding of Tanatar and Ceperley that the 2DEG be- 
comes a Wigner solid for = 37 ± 5. 

Using the Tq and the hard-disk bridge parameter 6*, 
we can calculate Ec{rsXiT) at any Our results (at 
T = 0) can be presented via a "polarization function" : 

^ i;.(r.,C)-gc(r.,0) ^ C^"^-d^°^-2 

a{r,) = Ci - C2lr, + C^/rl'^ - d/r]/^ (6) 

Here C± = (1±C)- In Hartree-Fock, the exponent a{rs) is 
a constant, a = {d + l)/d where d is the dimensionality. 
Hence, for the 2DEG, the Hartree-Fock value of a is 1.5. 
The coefiicients C1-C4 are 1.54039, 0.0305441, 0.296208, 
and 0.239047 respectively. A comparison of the 2-D and 
3-D a{rs) is given in Fig. 1(b). 

Our classical mapping of the 2DEG uses the GFMG 
spin-polarized xc-energy fit of Tanatar and Ceperley.El 
Our unpolarizcd energies (Table I) are in better agree- 
ment|-.with the more recent QMC results of Rapisarda 
et alE In Fig. 2 we compare our PDFs with those of 
GFMC at rs — 1, 5, and 20 given in Tanatar et al, and 
find satisfactory agreement. 

Equation (6) for p(rs,(^) pcpvides a comparison with 
the results of Varsano et al.p for the correlation ener- 
gies at partial polarizations, for 7's=20 and 30. Their 
-Ec{rs = 30, C) (an.) are, .01216, .01132, .00987, .00776 
at C=-2, .4, .6 and 0.8, while the results from our fit, eq. 
(6), are .01219, .01134, .00989, and .00775 respectively. 
Varsano et al. suggest that the fully spin-polarized phase 
is more stable, " by an amount close to the statistical un- 
certainty" . Instead of using the fit, we have evaluated the 
relevant PDFs accurately at rs=30 and T = 0. Varsano 
et al report a stabilization energy of -1.185 10~^ au., 
while ours is even smaller and not significant. 

When is large, the kinetic energy is less impor- 



tant; exchange favours spin-polarization and correlation 
favours the paramagnetic phase. Figure 3 displays the 
Hclmholtz free energy for rs= 5, 10, 20 and 30. Stable 
spin-polawzed phases in the 3DEG have been reported 
reccntlyJl3 Studies related to the metal-insulator transi- 
tion in the 2DEG, and numerical studies of 2-D clusters, 
have increasingly suggestecLthe possibility of ferromag- 
netic phases of the 2DEG.EZl The 2DEG metal-insulator 
transition is influenced by the presence of impurities. 
Our results are for the pure, uniform system and there 
seems to be no stable ferromagnetic phase except near 
rs=30. For 2DEGs in GaAs, = 30 and T/Ef=0.1 cor- 
responds to 4.7 10* electrons cm~^ and 0.019 Kelvin. The 
more favourable m* of Si brings the densities close to the 
~ 10^° electrons cm~^ range. Details of the finite-T cal- 
culations, the effect of impurities on the spin-polarized 
phase etc., will be presented elsewhere. 

The computer codes used here may be remotelyj-ac- 
cessed by interested workers by visiting our websiteEj 

In conclusion, we present a classical calculation of the 
PDFs, correlation energies etc., of the 2DEG for any 
spin-polarization and temperature, and find good agree- 
ment with quantum-Monte-Carlo results at T = 0. The 
temperature mapping of the quantum effects obtained 
for 3-D and 2-D systems is found to be universal, to 
within the accuracy of the available quantum Monte- 
Carlo exchange-correlations energies. New results for 
spin-dependent correlation energies and spin polarized 
distribution functions are given. A comparison of the 
free energies of the polarized and unpolarizcd phases at 
finite-T is presented. 
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TABLE I. Relevant parameters and a comparison of the 
QMC correlation energies (au.) of Tanatar and Ceperley 
(TC), Rapisarda and Senatore (RS), with those from the clas- 
sical map-HNC (CHNC) for the paramagnetic 2DEG. 



T/Tf 



V 



-£c(TC) -£c(RS) -£c(CHNC) 



1 2.000 0.500 0.235 0.10850 - - 0.10838 

5 1.849 2.704 0.372 0.04775 0.0489 0.04910 

10 1.682 5.946 0.425 0.03034 0.0302 0.03023 

20 1.468 13.62 0.469 0.01758 0.0174 0.01739 

30 1.331 22.54 0.486 0.01251 0.0124 0.01245 



FIG. 1. (a) The quantum temperature for the 2DEG 
(dashes) and the 3DEG (solid line) as a function of the classi- 
cal coupling constant F — P/ts. (b) The dependence of the 
2DEG (dashes) and 3DEG (solid line) polarization functions 
at T = 0. 



FIG. 2. Here the HNC g{r) are compared with GFMC 
simulations of Tanatar and Ceperley. Solid lines: HNC, boxes: 
GFMC, Panels (a),(b) unpolarized case,((^ = 0), rs=l, and 
Ts — 20. Panels (c),(d) rs=5 spin-polarized (^ — 1) and 
unpolarized systems. 



FIG. 3. The total Helmholtz free energy of the unpolar- 
ized {(, — 0) and polarized {( = 1) phases of the 2DEG near 
T = and for = 5, 10, 20, and 30, as a function of T/Tp. 
Solid lines: C — 0; dashes with data points (squares): C = 1- 
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